Mon. Not. R. Astron. Soc. 000, ??-[TT] (2011) Printed 14 January 2013 (MN WJ^ style file v2.2) 



AGN feedback and iron enrichment in the powerful radio 
galaxy, 4C+55.16 



o 

<: 

00 
(N 

O: 

u 

o 

> 

in 
oo 
in 

(N 

rn 
O 



J. Hlavacek-Larrondo^*, A. C. Fabian^, J. S. Sanders^ and G. B. Taylor^ 

^Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CBS OH A 
^Department of Physics and Astronomy, University of New-Mexico, Albuquerque, NM 87131, 



Accepted 2011 April 26. Received 2011 April 22; in original form 2011 March 13 



USA 



ABSTRACT 

We present a detailed X-ray analysis of 4C+55.16, an unusual and interesting radio 
galaxy, located at the centre of a cool core cluster of galaxies. 4C+55.16 is X-ray bright 
(Lx(ciustcr) ^ 10^^ erg ), radio powerful, and shows clear signs of interaction with 
the surrounding intracluster medium. By combining deep Chandra (100 ks) with 1.4 
GHz VLA observations, we find evidence of multiple outbursts from the central AGN, 
providing enough energy to offset cooling of the ICM (PbubWcs = 6.7 x 10*'' erg s~^ ). 
Furthermore, 4C-f 55.16 has an unusual intracluster iron distribution showing a plume- 
like feature rich in Fe L emission that runs along one of the X-ray cavities. The 
excess of iron associated with the plume is around lO^M©. The metal abundances 
are consistent with being Solar-like, indicating that both SNIa and SNII contribute 
to the enrichment. The plume and southern cavity form a region of cool metal-rich 
gas, and at the edge of this region, there is a clear discontinuity in temperature (from 
kT ~ 2.5 keV to kT - 5.0 keV), metalhcity (from - 0.4Zq to - O.8Z0), and surface 
brightness distribution, consistent with it being caused by a cold front. However, 
we also suggest that this discontinuity could be caused by cool metal-rich gas being 
uplifted from the central AGN along one of its X-ray cavities. 
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1 INTRODUCTION 

AGN feedback plays a major role in quenching star forma- 
tion, enriching the surrounding medium with metals, and 
fuelling the supermassive black hole (SMBH) of the host 
galaxy. Major advancements in understanding the details of 
how AGN feedback operates have been possible through de- 
tailed studies of X-ray cavit i es (see a review on the to pic by 
iMcNamara fc NulsenI [20071 : iPeterson fc FabianI lioo^ l . The 
SMBH lying at the centre inflates these cavities, also known 
as bubbles, through jets of relativistic plasma. Bubbles 
therefore provide a direct measurement of the energy being 
injected by the SMBH into the surrounding medium. 

In cool core clusters of galaxies, the active galactic 
nuclei (AGN) lying at the centre can energetically offset 
cooling of the intracluster medium by inflating these large 
cavities correlated with radio lobes, inducing weak shocks 
and propaga,ting energy through sound/pressure waves 
( Birzan et all |2004 iRaffertv et al] I2OO6I: [Punn fc FabianI 



This paper aims to study the AGN feedback processes 
arising in 4C-I-55.16, an unusual and interesting radio galaxy, 
located at the centre of a cool core cluster of galaj c ies at 
a r edshift of z = 0.2412 (iPearson fc Readheadlll98ll. Il988l: 
Wh yborn et~al]|l985l : iHutchings et al.lll988l : ISchneider et all 
i20(m . 4C-H55.16 is X-ray bright (Lx ~ 10*^ erg s'^ ), ra- 
dio powerful (I/R = 8Jy beam~^ at 1.4GHz), and shows 
clear signs of interaction between i ts central galajcy and 
the intracluster medium (IGM) (see IWhvborn et al.l Il985l : 
llwasawa et al. I ll999l . l200ll ). providing an interesting case for 
studying AGN feedback in clusters of galaxies. 



Furthermore, 4C+55.16 has an unusual i ntracluster iron 

i, llwasawa 



distribution. Using a 10 ks Chandra exposure,! 



2006, 2008': 'McNamara fc Nulscn"2007; Fabia n et all 12001 
2006. : .Forman et al.,,2005. : , Sanders fc Fabiani,2d07l ). 
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~et al.l 

(2001) found that there was a large increase in metallicity at 
a radius of about 10 arcsec (~ 40 kpc), that went from half 
solar to twice solar. The authors suggested that this increase 
was due to a plume-like structure located in the south-west 
side of the cluster, which had a strong Fe L emission feature. 
The reason as to how so much iron could be accumulated re- 
mained unclear. Here, we use significantly deeper Chandra 
observations of the source (100 ks) to study the abundance, 
and its distribution. By combining the Chandra observa- 
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Figure 1. Left - Exposure-corrected 0.5 — 7 keV Chandra image of 4C+55.16, spatially smoothed with a 1" gaussian function and 
covering ~ 500 X 500 kpc (or 130 X 130"). Shown in this image are the inner cavities, which are filled with 1.4 GHz radio emitting 
particles, the other possible cavity (with no GHz radio emission), the outer and inner edge as seen i n the X-ray , and the plume-like 
structure associated with the abundance increase. Right - Hubble Legacy Archive image of the cluster l|Edgell200oii . 



tions with VLA radio data, we make a detailed study of the 
X-ray cavities of the source. 

Section 2 presents the details of the data reduction. In 
Section 3 and Section 4, the techniques and results concern- 
ing the imaging and spectral analysis are respectively shown. 
Finally, we discuss the results in Section 5 and present the 
conclusions in Section 6. We adopt = 70 km s"'^ Mpc~^ 
with f2m = 0.3, r^A = 0.7 throughout this paper, and all er- 
ror estimates are on the Icr level. The redshift of the source 
corresponds to a scale of 3.807 kpc arcsec"^. 



2 DATA REDUCTION 
2.1 Chandra - X-ray 

4C-I-55.16 was originally observed with Chandra on 2000 
October 10 for 10 ks (ObsID 1645) with the Chandra CCD 
Imaging Spectrometer (ACIS) in FAINT mode, such that 
the cluster was centred on the ACIS-S3 back-illuminated 
chip. It was subsequently observed on 2004 January 3 for 100 
ks (ObsID 4940) in VFAINT mode, significantly improving 
the image quality, and was centred on the ACIS-S3 back- 
illuminated chip, with the ACIS-Sl, ACIS-S2, ACIS-I2 and 
ACIS-I3 also switched on. 

Both ObsID were processed, cleaned and calibrated 
using the latest version of the CIAO software (ciAOv4.3, 
CALDB4.4.1), and starting from the level 1 event file. We ap- 
plied both CTI (charge time interval) and time-dependent 
gain corrections, as well as removed fiares using the other 
back-illuminated chip (ACIS-Sl) and the lc_CLEAN script, 
with a 3(T threshold. 

For ObsID 4940, some strong background flares were 
seen towards the end of the exposure, and were removed. 
This resulted in a final exposure time of 73.6 ks. We then 
exposure-corrected the image, using an exposure map gener- 
ated with a monoenergetic distribution of source photons at 



1.5 keV(which is almost the peak energy of 4C-f 55.16). The 
exposure-corrected 0.5 — 7.0 keV image for ObsID 4940 is 
shown on the left side of Fig.[T] We also reduced the previous 
observations (ObsID 1645), and although some fluctuations 
were seen in the data (about 20 perc ent from the mean in t he 
light curve), we kept all data as did llwasawa et al] l|200ll ). 

Spectra w ere analysed using xspec (vl2.6.0d,f, e.g. 
lArnaudI Il996l ) . The Galactic absorptions througho ut this 
paper were kept frozen at the lKalberla et al] (|2005l ) value, 
4.29 X 10^° cm"^ (log(iVH) = 20.63). Letting the absorption 
free to vary for the spectrum of the entire cluster did not 
i mprove signiflcantly t he flt . We also used the abundances 
of I Anders fc Grevessd (|l989l ) throughout this paper. 



2.2 VLA - Radio 

There exists a short exposure (~ 6 mins) of published VLA 
observations at 1.4 GHz and bandwidth of 50 MHz in con- 
figuration A (Xu et al. 1995). This configuration is needed 
in using the VLA to produce maps with similar spatial res- 
olution to Chandra. We retrieved the data from the NRAO 
archive and reduced it in AlPS (astronomical image PRO- 
CESSING system) using the standard procedures. We used 
3C286 as a flux calibrator, and the difmap softwa re to clean 
and self-calibrate the data (|Shepherd et al.ll 19941 ). 

The VLA observations reveal a strong core and twin 
lobes. The shorter of the two lobes points to the north-west, 
while the longer lobe is oriented almost directly to the south . 
4C-I-55.16 has been imag ed with VLBI (jPolatidis et al.ll 19951 : 
IXu et al.|[l995l ) and shows a core and one-sided jet pointing 
to the north-west. 
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Figure 2. Left: Hubble Legacy Archive image of the central galaxy. This image shows evidence of an optical jet associated with the 
central point source. Middle: Smoothed Chandra image of 4C+55.16 in the 0.5 — 7.0 keV energy band, where the CIAO sub-pixel scripts 
have been applied. Right: Greyscale is an unsharp-masked Chandra image where a 2D gaussian smoothed image (cr = 10 pixels) was 
subtracted from a less smoothed image ((7 = 2 pixels), and which reveals deviations in the original Chandra image. Also shown in light 
blue are the contours of VLA data at 1.4 GHz , configuration A (see Xu et al., 1995) with a beam size of 2.3" X 2.3". Contours are drawn 
between 0.003Jy beam~^ to 0.403 Jy beam~^ (5 per cent of the peak intensity), in a log scale. At least two X-ray cavities are seen: one 
to the south, and one to the north-west of the central point source. However, there seems to be another cavity to the west of the cluster. 



3 IMAGING ANALYSIS 

We use the deep Chandra observations (100 ks) to analyse 
the physical state of the hot ICM, unless mentioned other- 
wise. Combining this data set with the 10 ks observations 
would not have improved significantly the signal-to-noise. 

In Fig. [T] (right), we show th e Hubble Legacy Archive 
image of the cluster (Edg el bOOOl ), and in Fig. [2j we show 
in same-scale images, the optical image taken from the 
Hubble Legacy Archive, the Chandra X-ray and unsharp- 
masked images, as well as the 1.4 GHz VLA radio contours 
of the central regions of the cluster. The Chandra image 
was obtained after applying the CIAO sub-pixel scripts in 
order to obtain images with the best resolution possible. 
The unsharp-masked image was obtained by subtracting a 
strongly 2D gaussian smoothed image (cr — 10 pixels or 
a — 4.9") from a less smoothed image (a = 2 pixels), and 
reveals deviations in the original Chandra image. At least 
two X-ray cavities are seen: one to the south, and one to the 
north-west of the central point source. However, there seems 
to be another cavity to the west of the cluster. The contours 
show that two radio lobes fill the X-ray cavities, consistent 
with the idea that these cavities are created by the central 
AGN. However, deeper observations would be needed to de- 
termine whether there is radio emission associated with the 
western cavity. 

In the optical, we can see the central galaxy surrounded 
by a diffuse envelope, along with a bright point source that 
coincides with the X-ray and radio point sources. We also 
see a small jet in the optical, where the axis of the jet is 
aligned with the radio filled X-ray cavities. There could be 
a counterpart of the jet in the X-rays, but the resolution is 
not good enough to resolve it (even with the CIAO sub-pixel 
scripts that give the best resolution for Chandra, see middle 
panel of Fig. (2)). 

To further enhance the surface brightness fluctuations 
in the X-ray image, at each point, we subtract from the X- 
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Figure 3. Fractional difference between each pixel and the aver- 
age of an ellipse centred on the X-ray point source and passing 
through the pixel. When constructing the ellipses, we use an aver- 
age ratio between the major (a) and minor (6) axis of the ellipses 
of a/b ~ 1.6, and an average rotated angle of the ellipses with 
respect to the north of 6 ~ 35°, counter-rotated. The image has 
then been smoothed with a cr = 2 pixels Gaussian. Also shown in 
light blue are the contours of VLA data at 1.4 GHz. 

ray image the average value within an ellipse centred on the 
X-ray point source. We used an ellipse rather than a cir- 
cle, since the surface brightness contours better follow this 
shape. Using the surface brightness contours, we determine 
the average ratio between the major (a) and minor (fe) axis 
of the ellipses, as well as the average rotated angle of the 
contours with respect to the north {9, counter-clockwise). 
We find a/b ~ 1.6 and d ~ 35°. Using these values, at each 
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point in the X-ray image, we build an ellipse passing through 
the point and centred on the central X-ray point source, and 
then subtract the average value within the ellipse and cal- 
culate the fractional difference. We then smooth the image 
with a 2D gaussian function of a = 2 pixels. This gives the 
image shown in Fig. [S] Initially, we remove all point sources, 
except for the one associated with the central galaxy. To do 
this, we identify them with the CIAO program wavdetect, 
and replace them with the average value of a surrounding 
background. 

Fig. [3] shows interesting structures, including bright 
rims surrounding the southern cavity, and a tail of enhanced 
material starting out from the southern cavity and plume- 
like feature, and making its way from east to north. 

As shown in the unsharp-masked image of Fig. [1] (right), 
there are at least two X-ray cavities (one to the south, and 
one to the north-west). However, if the potential cavity to 
the west is real, we should see a disturbance in the surface 
brightness distribution. In the top panel of Fig. U we show 
different regions we selected along each of the cavities. For 
each region, we determined the average count-rate pixel"^, 
and then plotted its distribution for each cavity as a function 
of radius in the lower panel of Fig. [4] This figure shows decre- 
ments in X-ray emission for each of the cavities, although 
it is quite difficult to tell for the north-western cavity. This 
cavity is however filled with radio emission, and we therefore 
still consider it as a cavity. 



4 SPECTRAL ANALYSIS 
4.1 Central AGN 

Using the deep Chandra observations, we first examine the 
nucleus of 4C+55.I6. Fig. [5] shows a bright X-ray point 
source located at the centre of the cluster, which coincides 
with the radio point source and optical point source of the 
central galaxy. Although the X-ray point source is bright and 
point-like, we estimate that there is no significant pileup. 
Pileup occurs when two or more photons are detected as on e 
event (see for more details lDaviij'200ll: lRu8sell et ai]l20ld ). 
The significance of pileup can be examined by comparing the 
amount of good grades (grades 0,2,3,4,6) to the bad grades 
(grades 1,5,7) for our point-like nucleus. We find the ratio 
of bad/good to be ~ 0.03. Typically, pileup starts becoming 
problematic when the fraction of bad grades exceeds 10 per 
cent of the good grades. 

To estimate the nuclear luminosity, we analyse the 
0.5 — 7 keV spectrum of the inner r = 1" region, centred 
on the central point source. We use a surrounding annu- 
lus (r = 2 — 3") as a background, and bin the spectra 
with a minimum of 30 counts per bin. We fit an absorbed 
power-law model to the background-subtracted spectrum, 
and include phabs absorption to account for Galactic ab- 
sorption. Initially, we only consider Galactic absorption and 
keep it frozen to the lKalberla et al] (|2005l ) value. We use y^- 
statistics to find the best fitting model, and let the photon 
index and normalisation free to vary. We find a photon index 
ofF = 1.55 ±0.17, and a reduced of 1.3. Using the CFLUX 
model in xspec, we obtain a flux estimate for the nucleus 
corrected for absorption, and then convert it into a luminos- 
ity using the luminosity distance. We find a 2 — 10 keV unab- 
sorbed luminosity for the nucleus of (1.1±0.2) x 10*^ erg s~^ , 
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Figure 4. Top: Unsharp-masked image shown in Fig. 1, with 
contours of regions used to look at the surface brightness distri- 
bution along each cavity. Bottom: Surface brightness distribution 
as a function of radius along each cavity. Also shown, are the 
locations of each cavity (between the dot-dashed lines). 



and that the value does not depend significantly on the pho- 
ton index to within ±0.3. There is also no evidence suggest- 
ing that the nuclear luminosity has changed within the 3 
year difference separating the observing dates of the lOOks 
and lOks data. 

If the Galactic absorption is also free to vary, we find 
that the model is no longer well constrained, and the esti- 
mated value for th e absorption is not c onsistent with the 
values predicted bv lKalberla et al] |200^ . Adding an addi- 
tional absorption at the redshift of the source does not im- 
prove the fit, and converges towards a null value for the ad- 
ditional absorption . We also fit a single-temperature mekal 
(|Mewe et al.ll 19951) model to the nuclear spectrum, and find 
a reduced of 1.5. The parameters are not well constrained 
and the estimated temperature (kT = lOls" keV) is quite 
large compared to the surrounding gas. 

If the X-ray surface brightness profile of the cluster 
scaled with radius as a power-law, then it would predict 
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a count rate of 30 counts pixel"^ for the nucleus. However, 
the X-ray point source is about three times as bright (90 
counts pixel"^). Yet, its intrinsic luminosity remains ~ 100 
times fainter than the total X-ray luminosity of the cluster 
( Lx ~ 10*^ erg s~^ ; see Hlavacek-Larrondo & Fabian 2011 
for a discussion of this issue). As seen in Fig. [51 the opti- 
cal shows evidence of a jet associated with the central point 
source. Since we are not able to resolve the jet in the X-rays, 
the nuclear luminosity we derive might have a beamed con- 
tribution from the jet, making the intrinsic AGN luminosity 
even fainter. 

4.2 Temperature and abundance maps 

To study the temperature and abundance distribution across 
the cluster, we bin different regions together using a Contour 
Binning algorithm which follows the surface brightness vari- 
ations (see Sanders ,2006). This technique is better suited 
to analyse the temperature and abundance variations along 
radial structures, as opposed to azimuthal structures. 

We begin by analysing the large scale structures within 
900 X 900 kpc. We bin the regions so that they have a signal- 
to- noise ratio of 50 (2500 counts), and restrict the lengths to 
be at most twice the widths. We extract a 0.5 — 7 keV spec- 
trum for each region, and use a region further away, but still 
within the chip, as a background. We then fit an absorbed 
MEKAL model to the data, and use C-statistics. We let the 
temperature, absorption and normalisation parameters free 
to vary. The resulting temperature and metallicity maps are 
shown in Fig. [S] For the temperature, the errors on each 
value vary from 6 per cent in the inner regions to ~ 30 
per cent in the outer regions, and for the abundance they 
vary from ~ 20 per cent (inner) to ~ 40 per cent (outer) . 

The large-scale temperature map shows some structure, 
especially a temperature jump between the cluster core (in- 
ner 50 kpc) and surrounding gas on the southern side. The 
metallicity map also shows some structure. To the north- 
east, we see a first metal rich region (IZq) located at about 
100 kpc, and then a second (0.6^©) at about 160 kpc. Fur- 
ther away (r > 200 kpc), it seems that the north-east side 
is more metal poor than the south-western side. 

Next, we analyse the inner regions of the cluster, within 
190 X 190 kpc. We bin the regions so that they have a signal- 
to-noise ratio of 30 ( 900 counts), and also use C-statistics. 
The resulting temperature and abundance maps are shown 
in Fig. [S] Here, the errors on each value vary from ~ 9 per 
cent in the inner regions to ~ 20 per cent in the outer regions 
for the temperature, and from ~ 30 per cent (inner) to ^ 60 
per cent (outer) for the abundance. 

The temperature map now reveals more clearly the tem- 
perature jump seen between the core and surrounding re- 
gion just below the southern cavity and plume-like feature. 
The temperature jumps by a factor of two, from colder to 
warmer, and is similar to those seen in cold fronts. In the 
temperature map, we see a second jump in temperature to 
the west of the cluster at about r 120 kpc. This time, the 
temperature jumps from a colder (fcT ~ 4 keV) to warmer 
region {kT ~ 8 keV) (see also Section 4.3.2). The metallicity 
map also clearly reveals that the plume-like feature (located 
to the south-west at r ~ 1 — 10") is metal rich when com- 
pared to a region within the same radius but excluding the 
plume (see also Section 4.3.4). 
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Figure 5. Large scale temperature (top) and abundance (bot- 
tom) map of the cluster, as derived using the Contour Binning 
algorithm (see Section 4.2). We binned regions within 900 X 900 
kpc, with a signal-to-noise of 50. The temperature units are keV 
and those of the metaUicity are Zq. 

4.3 Spectral profiles 

4-3.1 Previous Chandra observations (lOks) 

Using the 10 ks Chandra observations. Ilwasawa et all (|200ll ) 
found that there was an large increase in metallicity at a 
radius of about 10 arcsec (~ 40 kpc), that went from half 
solar to twice solar, which they associated with the plume- 
like feature. 

In order to confirm this increase in metallicity, we re- 
duce the lOks observations using the latest Chandra reduc- 
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Figure 6. Temperature (left) and abundance (middle) maps of the inner regions, as derived using the Contour Binning algorithm (see 
Section 4.2), as well as exposure-corrected image of the inner regions of the cluster. We binned regions within 190 X 190 kpc, with a 
signal-to-noise of 30. The temperature units arc kcV and those of the mctallicity are Zq. 



tion scripts. We then proceed to a nalyse the data in the 
same way as in llwasawa et al] (12001 ). but choose to use the 
value derived by iKalberla et al. I (|2005l ) for the Galactic ab- 
sorption (4.29 x 10 '^ ° cm~ '^) instead of the one derived by 
iDickev fc LockmanI (|l990l . 4.20 x 10^° cm"^). Since these 
values are almost the same, we do not expect our results to 
var y significantly based o n this choice of parameter. 

Ilwasawa et all ll200ll ) areued that the redshift inferred 
from optical spectroscopy was not consistent with the one 
inferred by the Fe K line. From their analysis of the Fe K 
line, they derived a redshift of z ~ 0.254tQ gj!,g. However, our 
spectra indicates otherwise. Using the deep 100 ks data, we 
extract a 0.5 — 7 keV spectrum within a radius of r = 1 — 9", 
since this region is known to contain a strong Fe K line (see 
also Fig. [8]). We then fit an absorbed mekal model to the 
spectrum, and let the redshift, temperature, abundance and 
normalisation parameter free to vary. We use statistics, 
and find a redshift of 2 = 0.238 ± 0.007, consistent with the 
value inferred by optical spectroscopy, and therefore choose 
to use the optically derived redshift (2 = 0.2412) throughout 
this paper. 

We now proceed to take spectra from seven annuli, 
follow ing the spectral profiles extracted in Ilwasawa et al.l 
('2001') (r = 1-2.5,2.5-5,5-9,9-15,15-25,25-40,40- 
60" j^. These annuli exclude the nucleus. We then fit an ab- 
sorbed MEKAL model to each spectrum, and derive the tem- 
perature, abundance and normalisation parameter using 
statistics. The results are shown in the left panels of Fig. 
[T] We do the same using the deep Chandra 100 ks obser- 
vations, but since the data is much deeper, we also include 
deprojection effects. The right panels of Fig. [7| show the re- 
sults we obtain for both the projected and deprojected data. 
The latter are obt ained following the Direct Spectral Depro- 
jection method of ISanders fc FabianI (|2007l ). 



As shown in Fig. [T] the 100 ks observations give much 
tighter constraints on the temperature and abundance pro- 
files, and most importantly they do not show the abundance 
increase seen in the 10 ks observations. The abundance in- 
crease we derive from the 10 ks observations is also not the 
same as in llwasawa et al.1 (|200ll ). who found an increase that 
went from half solar to twice solar. We only find an increase 
form half solar to one solar. However, when analysing in 
more detail the 10 ks observations, we first notice that the 
model is not very sensitive to the abundance. Even if we 
freeze the abundance to four times solar, it only increases 
the from 40.15 to 42.19. The abundance increase we see 
from the 10 ks observations should therefore be taken lightly. 

Secondly, wc analyse the region between r = 1 — 9", 
as did Iwasawa ct al. (20(}l|), but using the deep 100 ks ob- 
servations. We fit an absorbed mekal model to the spec- 
trum, and derive a best fitting value of kT = 2.93 ±0.11 
keV and Z = 0.74 ± 0.12 solar. When analysing this re- 
gion, Ilwasawa et al.l (120011 ) found using the 10 ks observa- 
tions that kT = 3.34 ± 0.35 keV and Z = l-^Sto ll solar. 
In Fig|8] we plot our spectrum, as derived from the 100 ks 
observations. Overlaid on top, we plot our best fitting model 
(black curve), and a mo del where we kept kT and Z frozen 
to the values derived by 'iwa sawa et ah! (|200ll ). and only al- 
lowed the normalisation parameter fre e to vary. As shown in 
the figure, the parameters derived bv Ilwasawa et ah! (|200ll ) 
seem to overestimate the strength of the Fe K line, and 
also underestimate the soft energy part of the spectrum. 
We suspect the latter may be caused by the new calibra- 
tion for the absorption mod el of the detector (see for more 
details iMarshall et al.ll2004l ). which now better corrects for 
the ACIS Chandra contaminant at soft energies. Hence, it 
is possible that the earlier reduction scripts over corrected 
the spectrum at soft X-rays, and allowed the fit to converge 
toward higher abundances. 



Ilwasawa et al. 1 1I2OOII I used eight profiles instead of seven, split- 
ting the first bin into two bins of r = 1 — 2.5" and r = 2.5 — 5". 
However, our fit could not converge for each of these bins sepa- 
rately, and we were forced to merge them together in order to get 
a good enough S/N ratio and derive parameters. 



4-3.2 Deep Chandra observations (lOOks) 

We now proceed to do a more in depth analysis of the 100 k s 
observations. The regions selected bv Ilwasawa et ah! (|200ll ) 
each had a different number of counts, ranging from ~ 2000 
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Figure 7. Temperature (top) an d metallicity ( b otton i) pro- 
files for tlie regions as selected in llwasawa et al.l 1 I2OOII ) (r = 
1 - 2.5, 2.5 - 5, 5 - 9, 9 - 15, 15 - 25, 25 - 40, 40 - 60"). The left 
panels show the results for the 10 ks observations, reduced us- 
ing the latest Chandra calibrations, and those on the right show 
the results derived from the 100 ks observations (projected and 
deprojected). The deeper data do not show the inc rease in metal- 
licity seen in the 10 ks observations and reported in llwasawa et all 
1I2OOII '). 




Energy (keV) 

Figure 8. Spectrum of the central r = 1 — 9" region (excluding 
the nucleus). The black curve shows the best-fitting model ob- 
tained for an absorbed (Galactic) MEKAL model (fcT = 2.93±0.11 
keV, Z = 0.74 it 0.12 solar). The red curve shows the best-fitting 
model while keeping the temperatu r e and abundance frozen to 
the values found by llwasawa et al.l l|200ll ) for the same region 
{kT = 3.34 ± 0.35 keV, Z = 1.93lg;g2 solar). 
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Figure 9. Projected (filled symbols) and deprojected (non-filled 
symbols) temperature and metallicity profiles, as well as depro- 
jected electron density, electron pressure and entropy profiles of 
the cluster, as derived by selecting annuli with a minimal signal- 
to- noise ratio of 63 (~ 4000 counts). 



counts to ~ 10000 counts. Instead, we choose to select re- 
gions containing roughly the same signal-to-noise (S /N~ 63 
or ~ 4000 counts) , and push the analysis to larger radii. We 
selected a background region within the same chip (ACIS- 
S3) but far from the cluster. The net count rate per pixel^ 
of the selected background region is also the same as the 
count rate per pixel^ of the ACIS-Sl chip (which is also 
back- illuminated), and should therefore not contain much 
cluster emission. 

The projected and deprojected temperature and metal- 
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licity, as well as deprojected electron density, electron pres- 
sure and entropy profiles {S = kTrie^^^) as a function of 
radius are show in Fig. [9] For each spectrum, we fitted an 
absorbed mekal model, and used statistics. Although we 
include in Fig.|9]the deprojected results, we found that they 
were consistent with the projected ones to within la error 
estimates. Two-temperature models were also not needed, 
as they did not improve the fit. 

Fig. O shows interesting features at r > 100 kpc. The 
temperature structure remains constant at large radii, at 
least within the error measurements, while the metallicity 
structure on average decreases with radius (from Z ^ 0.75 
solar to Z 0.2 solar). 

There seems to be a slight increase in metallicity at 
r ~ 160 kpc. This could be due to the same feature no- 
ticed earlier while analysing the large-scale metallicity map 
(see bottom panel of Fig. [5]), which showed a slight increase 
at a similar radius on the north-eastern side of the cluster. 
However, the increase in metallicity is well within the Icr un- 
certainties, and therefore does not appear to be statistically 
significant. If we compare the spectra of the different radial 
intervals with one another, we find no significant change in 
the iron line complexes. 

However, there appears to be an increase in temper- 
ature at r ~ 130 kpc (Fig. |9]), accompanied by a break in 
electron density, electron pressure and entropy, to more than 
a Icr level. This could be associated with the warmer region 
noticed earlier in the temperature map of the inner regions 
(Fig. 6), where we could see a jump in temperature from 
colder [kT ~ 4 keV) to warmer \kT ~ 8 keV). 

4-3.3 southern X-ray cavity 

We now concentrate on the spectral analysis along the south- 
ern cavity. Clear bright rims can be seen surrounding this 
cavity, and could indicate the presence of a shock front (see 
Fig. [3]). Using regions along the southern cavity, we extract 
a spectrum, and fit an absorbed mekal model to each re- 
gion. The absorption is again taken as the purely Galactic. 
The resulting temperature profile is shown in Fig. IIQI This 
figure shows that there is a temperature jump associated 
with the edge of the cavity at r ~ 40 kpc, but it does not 
indicate a shock front. Instead it shows a temperature jump 
going from a cooler to warmer region (from kT ^ 2.8 keV 
to kT ~ 5.0 keV), i.e. it could be a cold front. In Section 
5.2, we discuss in greater detail this front. 



Figure 10. Projected temperature profiles along the southern 
cavity (red) and plume-like feature (black) , derived using a single 
MEKAL model, as well as electron density, electron pressure and 
entropy profiles of the plume-like feature. The potential cold front 
is located at r ~ 40 kpc. Selected regions are shown in the top 
figure which is a 3 colour image (red 0.3 — 0.7, green 0.7 — 1.0, 
and blue 1.0 — 7.0 keV). The metal-rich plume stands out in green 
as it is dominated by the Fe L emission feature, centred at ~ 0.8 
keV. 



4-3.4 Plume-like feature 

According to the middle panel of Fig. [6l there is an increase 
in metallicity associated with the plume-like feature. We can 
also see from the left panel of Fig.[6l that there is an increase 
of temperature associated with the edge of the plume-like 
feature. To be more precise, in Fig. 1101 we plot the temper- 
ature profile, electron density and pressure, as well as the 
count rate per pixel'^ (equivalent to the surface brightness 
distribution) as derived by selecting regions along the plume- 
like feature. However, we derive a count rate per pixel'^ for 
every arcsec in order to better trace the potential distur- 
bances in the surface brightness profile. This figure shows 
a temperature jump at a radius similar to the one associ- 
ated with the southern cavity (r ~ 40 kpc). We associate 
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Figure 11. Ratio of the spectrum from the region containing 
the plume-like feature (filled points) to the best-fitting MEKAL 
model of the spectrum not containing the plume feature (non- 
filled points), both within r = 1 — 9". The normalization of the 
region containing the plume was determined by fitting the spec- 
trum excluding data between 0.7— 1 keV. This plot clearly shows 
the excess of the Fe L emission feature associated with the plume. 



this temperature jump with the same cold front that can be 
seen at the edge of the southern cavity, and which we discuss 
in Section 5.2. 

We now compare the spectrum of the plume to the one 
of a region of same radius {r — 1 — 9"), but excluding the 
plume. The region containing the plume is taken as being 
between 290° and 340° (counter-clockwise from the north), 
and the remaining region is taken as the one excluding the 
plume (i.e. a C-shape). In Fig. llll we plot the ratio between 
the spectrum of the region containing the plume- like feature 
(filled points) to the best-fitting mekal model of the spec- 
trum not containing the plume feature (non-filled points). 
We used statistics and the normalization of the region 
containing the plume is determined by fitting the spectrum 
excluding data between 0.7—1 keV. The excess of the Fe L 
emission feature associated with the plume is evident. 

This excess could nonetheless be caused by the colder 
temperature material associated with the cold front (and 
which coincides with the plume- like feature) . However, even 
if we let the temperature free to vary, the ratio between 
the spectrum of the region containing the plume-like fea- 
ture (filled points) to the best-fitting mekal model of the 
spectrum not containing the plume feature still shows an 
excess in the Fe L emission feature. The excess is therefore 
still likely to be caused by a metallicity effect. 

Its nature could b e due t o iron enrichment by SNe, as 
suggested by Iwasaw a et al.l l|200ll ). We can test this hy- 
pothesis by using the variable abundance model vmekal in 
XSPEC. Using the deep 100 ks Chandra data, we analyse 
the spectrum of the plume-like feature within r = 1 — 9" 
(see Fig. |8|. The best-fitting mekal model finds a temper- 
ature of kT = 2.93 ± 0.11 keV and metal enrichment of 
Z = 0.74 ± 0.12 solar (equivalent to Zpe as the spectrum is 
dominated by iron emission lines). We also find a reduced 
statistic of 1.26. 



If we test the vmekal model using the enrichment 
pattern by SNIa, follow ing the abundance ratios given by 
iDupke fc ArnaudI (|200ll ). we find kT = 2.99 ±0.11 keV, 
Zfc ~ 0.47 ± 0.07 solar and a reduced statistic of 1.37 
(same number of degrees of freedom). On the other hand, 
using the enrichment pat tern by SNII, also follo wing the 
abundance ratios given bv lDupke fc ArnaudI (|200ll ). we find 
kT = 2.90 ± 0.14 keV, Zpe = 0.25 ± 0.06 solar and a re- 
duced statistic of 2.02. Hence, the enrichment pattern by 
SNIa provides a much better fit, but not quite as much as 
the Solar photosphere enrichment pattern. We discuss these 
results in Section 5.3. 



5 DISCUSSION 
5.1 X-ray cavities 

There are at least two, and possibly three, X-ray cavities 
associated with the central AGN. The energy stored within 
each of the cavities can be es timated u sing Eg. [1 (e.g . 
Birzan et al] 120041: iDunn et all [2005: Raffertv et alTl2006l : 
bunn fc Fabianll2006l . |2008| ). Here, P is the thermal pres- 
sure of the ICM at the radius of the bubble and estimated 
from the X-ray data, V is the volume of the cavity and for 
a relativistic fiuid 71 = 4/3, therefore -Ebubbio = 4PV (this 
is also supported observationally, see [Graham et al.ii200S ). 



-Ebub 



71 



-PV 



71 



(1) 



We assume that the cavities are of ellipsoidal shape. 
The volume is then given hy V = 47vRl,Ri/3, where Ri is 
the semi-major axis along the direction of the jet, and i?w 
is the semi-major axis perpendicular to the direction of the 
jet. 

The power injected into the medium is determined 
by dividing the energy of the bubble with its age. 
The latter is given by the buoyant rise time, the refill 
time or the sound cr ossing t ime. See respectively E q. [2] 
JChurazov et all 1200 ll). Eg. [51 (iMcNamara et al.ll2000l ) and 



Eg. [41 (|McNamara fc Nulsenfl2007l )T Here, R is the distance 
from the radio point source to the middle of the cav- 
ity (projected), S is the cross-sectional are a of the bubble 
(S = ttJ?^), Cd=0.75 is the drag coefficient jChurazov et al.l 
l200ll ). g is the local gravitational acceleration [g = GM{< 
R)/R^), r is the bubble radius (r — (RiR^)^^'^ for an el- 
lipsoidal bubble), and Cg is the sound crossing time (cs = 
^y 72fcr/ ( /^m jj), where kT is the plasma temperature at 
the radius of the bubble, 72 = 5/3 and fi — 0.62). 



uoyant 



R 



SCp 
2gV 



t 

''CS 

Cs 



(2) 



(3) 



(4) 



The buoyant rise time is the time it takes a bubble to 
reach its terminal buoyant velocity, which depends on the 
medium drag forces. This is a good estimate of a bubble's 
age that has clearly detached from their AGN and has risen. 
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such as the cavity to the west in 4C+55.16. The refill time 
is the time it takes a bubble to rise buoyantly through its 
own diameter starting from rest (which is probably not the 
case). The sound crossing time is the time it takes a bubble 
to travel at the speed of sound. The latter is used under 
the assumption that the bubbles travel at subsonic speeds. 
Between the three time estimates, it is still not clear which is 
the best to use, although they should not vary significantly. 

Table[T]shows the estimates we obtain for the power be- 
ing injected into the ICM by each of the bubbles. 4C+55.16 
has a cooling luminosity in the 0.5 — 7 keV range of Lcooi = 
(1.99 ± 0.02) X lO"** erg s"^ (rcooi = 45 kpc for a cooling 
time of 3 Gyr). The power being injected into the medium 
by all three cavities is therefore sufficient to prevent the gas 
fro m cooling, at l e ast w ithin a factor of two. T he studies 
by iRaffertv etahl (120061 ) and iDunn et al.1 l|2005h included 
4C-I-55.16 in their sample of clusters with X-ray cavities. 
Our results are consistent with theirs, on average within a 
factor of two. 



5.2 Cold front 

Cold fronts in cool core clusters are contact discontinuities 
thought to originate from sloshing of low-entropy gas in the 
cluster core caused by disturbances on the central poten- 
tial by ph enomena such as past subclus ter mergers (see a 
review by iMarkevitch fc Vikhlininll2007l ). Simulations sup- 
port this idea, but require that clusters with cold fronts 
have steep entropy profiles (e.g. lAscasibar fc MarkevitdJ 
I2006D . This is generally the case for cool core clusters. Al- 
though cold fronts could harbour significant amounts of ki- 
netic ene rgy, it is not clear that they could dissipate it ef- 
ficiently (|Markevitch et all 1200 ll ). The exact cause of the 
sloshing is also not well known. Past mergers could cause 
slos hing of the cluster core (as supported by simulations, 
see lAscasibar fc Markevitch 2006), but another possibility 
is whether outbursts from the central AGN could cause the 
core to recoil and slosh. 

In Section 4., we showed that 4C-I-55.16 has a clear tem- 
perature jump on the south-western side of the cluster asso- 
ciated with the southern cavity and plume-like feature, and 
that the plume-like feature is about twice as rich in metals 
than a region located within the same radius but excluding 
the plume (see the temperature and abundance map in Fig. 
[S| . A more detailed analysis also showed that there is a dis- 
continuity in the surface brightness profile associated with 
this front, see bottom panel of Fig. 1101 (although the break is 
at a radius slightly further away from the jump in tempera- 
ture). These results are all consistent with the discontinuity 
being a cold front. 

In the literature, we find several examples of cold 
fronts showing a metallicity jump like in 4C-f55.16, but 
with no cl ear dis continuit y in the temperature profile, 
e.g. A2052 llde Pla,a ct al . 201oh . A2199 and 2A 0335+096 
l|Sanders fc Fabianll2006bl ). Others do not show any evidence 
of a metallicity jump, but clearly show a discon tinuity in 
temperature, such as A496 (|Dupke fc Whitell2003l ') or A2204 
l|Sanders et al.ll2005l ). Our front also shows a break in en- 
tropy and pressure. The pressure remains constant across 
the front, a typical feature of cold fronts. Another key point 
that seems to support the idea that the front we observe 
is a cold front is the spiral (or tail-like) structure we see 



in Fig. O This figure was obtained by subtracting at every 
pixel, the average value of an ellipse passing through the 
point, centred on the X-ray point source. The figure shows 
an excess of emission (~ 30 per cent from the mean) in the 
shape of a large one-armed spiral structure extending from 
the metal-rich plume feature to the outskirts of the cluster. 
The mass within this structure is roughly ~ 6 x 10^^ Mq, if 
we assume a cylindrical shape. 

Similar spiral-like structu res have been seen in other 
clusters, s uch as in A2029 (IClarke et ahl 120041 ') . in 2 A 
0335-1-096 (ISanders et al. 20091) and ev en in Perseus (e.g. 
IChurazov et al.ll2003l : lFabian et al.ll2000l ). Simulations in the 
context of cold fronts have shown that gaseous merging 
clusters with non zero impact parameters are able to re- 
produce these spiral structures as the central cool gas ac- 
quires angular momentum, in addition to the cold front (e.g. 
lAscasibar fc Markevitchll2006l ) . It takes several Gyrs for the 
cluster to regain a relaxed morphology. We expect 4C-I-55.16 
to be relaxed, since the cluster shows no obvious signs of re- 
cent merger activity. The X-ray cavity to the south could 
therefore still have has the time to form, post merger, as 
bubble time scales are much shorter (~ 0.1 Gyr). 

If the discontinuities we see in temperature, abundance, 
pressure and entropy are not caused by a merger-induced 
cold front, they could simply by caused by feedback of the 
central AGN. As bubbles rise, the y are expected to dra g 
behind them, cool, metal-rich gas (|Churazov et al.ll200ll ). 
However, if there is originally cool low-entropy metal-rich 
gas at the base of the bubble, as its being inflated by the 
AGN, it could push the cooler metal-rich gas outwards. The 
gas would eventually slide back down around the bubble 
and would cause the appearance of a plume-like metal-rich 
feature running along the cavity, just like in 4C-I-55.16. 

5.3 Metal enrichment 

We now focus on analysing the metal enrichment associ- 
ated with the plume-like structure. This structure is rich 
in the Fe emission feature (see Fig. Ilip . It is about QAZq 
more metal rich than a region located within the same ra- 
dius, but excluding the plume (which has a metallicity of 
~ 0.45Zq). If we assume that the plume has a cylindrical 
shape with a length of 30 kpc and width of 14 kpc, and that 
the density is ~ 0.065 cm~^ (based are the deprojected value 
along the plume), then the excess of iron mass is around 
8.9 X 10^ Mq. It is worth mentioning that our results are 
based on the fits we have obtaine d using the abundance ra- 
tios of lAnders fc Grevesse' ('l989l). whereas the more recent 
tables of Asplund ct al. (200i) or iLoddersI (|2003l ) predict 
an iron abundance about 1.5 times lower. Refitting the data 
with these tables instead, we find that the plume has an iron 
excess of ~ 0.6^© instead of ~ OAZq. However, the excess 
of iron in terms of mass does not change, as this quantity 
also depends on the relative abundance ratios, and the net 
effect cancels out. 

We can estimate the age of this structure in several 
ways. First, by using a measure of th e diffusion coefficient , 
such as the one derived for Perseus in lRebusco et al" I (l2005l . 
2 X 10^^ cm'^s"^), and an average size of the plume (30 x 14 
kpc^), we find that time-scale associated with the plume 
should be on the order of ~ 6 x 10^ yr. Second, we can 
estimate the age of the plume if we consider that it has to 
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9.9 
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0.80 
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be younger than the local cooling time. The temperature 
and density of the plume then correspond to a cooling time 
of ~ 5 X 10* yr. We are ignoring heating here, which may 
be plausible since turbulence has not yet had the time to 
destroy it according to the long diffusion time. An age of 
about 5 — 6 X 10* yr implies that if the iron enrichment of 
the plume was caused by SNIa, then it would require around 
2.5 SNIa per century, if we assume that 0.7Mq are produced 
for every SNIa. This is at least an order of magnitude higher 
than the predicted rate of SNIa in an elli ptical galaxy of 
similar Lb (see e.g. ICappellaro et al.lll997l ). 

The ages we derive for the plume are in agreement with 
those of the X-ray cavities (tbubbic ~ 10^ - 10* yr, see Table 
1), and are consistent with the idea that the plume could 
be created by the uprising of the southern bubble. Other 
examples of clusters have also shown evidence of significant 
quantities of metal-rich gas (10* — IO^Mq) being uplifted 
by buoyantly rising bubbl es. These include S ersic 159-03 
JWerner et al.ll201lli, M87 l| Werner et al.ll2010l ) and Hydra 
A i Simionescu et al]|2009l ). In order for the plume- like struc- 
ture to be long lived, magnetic fields could contribute to sta- 
bilizing it s uch as in the case o f the emission-line fila ments in 
Cent aurus l|Tavlor et al.ll2007l ) and in Perseus (Fabi ajTet al.l 
l2008l ). 

In Section 4, we analysed the enrichment pattern of the 
plume using either SNIa or SNII enrichment s. Past studies, 
including the previous one on 4C-I-55.16 bv llwasawa e t al] 
l|200ll ). suggested that the central ir on excess seen in c luster s 
was only caused by SNIa (see also [Matsushita et al]|2003h . 
However, our results led us to conclude that although the 
spectrum is more consistent with it being enriched by SNIa 
compared to SNII, the solar enrichment pattern still pro- 
vides a better fit. Sa nders fc Fab ian (2006a) also find that 
the enrichment pattern of the core of the Centaurus cluster 
is more consistent with Solar values, suggesting that both 
SNIa (~ 70 per cent) and SNII (~ 30 per cent) contribute to 
the e nrichment (see also Ide Plaa et al.ll2007l : iLovisari et al.l 
Our results agree with this picture. 



6 CONCLUDING REMARKS 

4C-I-55.16 is X-ray bright ( Lx ~ lO**^ erg s~^ ), radio pow- 
erful, and shows clear signs of interaction between its central 
galaxy and the ICM. 4C-I-55.16 has at least two, and possi- 
bly three X-ray cavities, two of which are filled with radio 
emitting particles. The power stored within them is around 
6.7 X 10'*'' erg s~^ , and is sufficient to prevent the ICM from 
cooling. Our study confirms earlier results suggesting that 
there is a plume-like feature, running along one of the cavi- 
ties, rich in the Fe L emission feature. However, we also find 
that the plume and cavity form a region of cool metal-rich 



gas that has a temperature and metallicity jump (by a factor 
of 2), and density jump. It could therefore be a cold front, 
or this could also be an example of metal enrichment by the 
central AGN, which has the potential to uplift cool metal- 
rich gas from the central galaxy. Finally, we find that the 
plume has an excess of iron (Mpo = 8.9 x 10^ Mq), and that 
its enrichment is more consistent with being Solar-like, and 
therefore suggesting that both SNIa and SNII contribute. 
This is in disagreement with the more general view that 
only SNIa enrich in the central regions of clusters. 



ACKNOWLEDGMENTS 

JHL recognizes all the support given by the Cambridge 
Trusts, Natural Sciences and Engineering Research Coun- 
cil of Canada (NSERC), as well as the Fonds Quebecois de 
la Recherche sur la Nature et les Technologies (FQRNT). 
ACF thanks the Royal Society for support. GBT acknowl- 
edges support for this work from the National Aeronautics 
and Space Administration through Chandra Award Num- 
ber GO0-11139X. 



REFERENCES 

Anders E., Grevesse N., 1989, Geochim. Cosmochim. Acta, 
53, 197 

Arnaud K. A., 1996, in G. H. Jacoby & J. Barnes ed.. Astro- 
nomical Data Analysis Software and Systems V Vol. 101 
of Astronomical Society of the Pacific Conference Series, 
XSPEC: The First Ten Years, pp 17-+ 

Ascasibar Y., Markevitch M., 2006, ApJ, 650, 102 

Asplund M., Grevesse N., Sauval A. J., Scott P., 2009, 
ARA&A, 47, 481 

Birzan L., Rafferty D. A., McNamara B. R., Wise M. W., 
Nulsen P. E. J., 2004, ApJ, 607, 800 

Cappellaro E., Turatto M., Tsvetkov D. Y., Bartunov O. S., 
PoUas C, Evans R., Hamuy M., 1997, A&A, 322, 431 

Churazov E., Briiggen M., Kaiser C. R., Bohringer H., For- 
man W., 2001, ApJ, 554, 261 

Churazov E., Forman W., Jones C, Bohringer H., 2003, 
ApJ, 590, 225 

Clarke T. E., Blanton E. L., Sarazin C. L., 2004, ApJ, 616, 
178 

Davis J. E., 2001, ApJ, 562, 575 

de Plaa J., Werner N., Bleeker J. A. M., Vink J., Kaastra 

J. S., Mendez M., 2007, A&A, 465, 345 
de Plaa J., Werner N., Simionescu A., Kaastra J. S., Grange 

Y. G., Vink J., 2010, A&A, 523, A81+ 
Dickey J. M., Lockman F. J., 1990, ARA&A, 28, 215 
Dunn R. J. H., Fabian A. C, 2006, MNRAS, 373, 959 
Dunn R. J. H., Fabian A. C, 2008, MNRAS, 385, 757 



© 2011 RAS, MNRAS 000, ??-[II] 



12 J. Hlavacek-Larrondo, et al. 



Dunn R. J. H., Fabian A. C, Taylor G. B., 2005, MNRAS, 
364, 1343 

Dupke R, White III R. E., 2003, ApJ, 583, L13 

Dupke R. A., Arnaud K. A., 2001, ApJ, 548, 141 

Edge A., 2000, in HST Proposal A Continuation of a 

SNAPSHOT survey of X- ray selected central cluster 

galaxies, pp 8719 — h 
Fabian A. C, Johnstone R. M., Sanders J. S., Conselice 

C. J., Crawford C. S., Gallagher III J. S., Zweibel E., 

2008, Nat, 454, 968 
Fabian A. C, Sanders J. S., Allen S. W., Crawford C. S., 

Iwasawa K., Johnstone R. M., Schmidt R. W., Taylor 

G. B., 2003, MNRAS, 344, L43 
Fabian A. C, Sanders J. S., Ettori S., Taylor G. B., Allen 

S. W., Crawford C. S., Iwasawa K., Johnstone R. M., Ogle 

P. M., 2000, MNRAS, 318, L65 
Fabian A. C, Sanders J. S., Taylor G. B., Allen S. W., 

Crawford C. S., Johnstone R. M., Iwasawa K., 2006, MN- 
RAS, 366, 417 

Forman W., Nulsen P., Heinz S., Owen F., Eilek J., 
Vikhhnin A., Markevitch M., Kraft R., Churazov E., 
Jones C, 2005, ApJ, 635, 894 

Graham J., Fabian A. C, Sanders J. S., 2008, MNRAS, 
386, 278 

Hlavacek-Larrondo J., Fabian A. C, 2011, MNRAS, pp 
138— H 

Hutchings J. B., Johnson I., Pyke R., 1988, ApJS, 66, 361 
Iwasawa K., Allen S. W., Fabian A. C, Edge A. C, Ettori 

S., 1999, MNRAS, 306, 467 
Iwasawa K., Fabian A. C, Allen S. W., Ettori S., 2001, 

MNRAS, 328, L5 
Kalberla P. M. W., Burton W. B., Hartmann D., Arnal 

E. M., Bajaja E., Morras R., Poppel W. G. L., 2005, A&A, 

440, 775 

Lodders K., 2003, ApJ, 591, 1220 

Lovisari L., Schindler S., Kapferer W., 2011, ArXiv e-prints 
Markevitch M., Vikhlinin A., 2007, Phys. Rep., 443, 1 
Markevitch M., Vikhhnin A., Mazzotta P., 2001, ApJ, 562, 
L153 

Marshall H. L., Tennant A., Grant C. E., Hitchcock A. P., 
O'Dell S. L., Plucinsky P. P., 2004, in K. A. Flanagan & 
O. H. W. Siegmund ed.. Society of Photo-Optical Instru- 
mentation Engineers (SPIE) Conference Series Vol. 5165 
of Presented at the Society of Photo-Optical Instrumen- 
tation Engineers (SPIE) Conference, Composition of the 
Chandra ACIS contaminant, pp 497-508 

Matsushita K., Finoguenov A., Bohringer H., 2003, A&A, 
401, 443 

McNamara B. R., Nulsen P. E. J., 2007, ARA&A, 45, 117 
McNamara B. R., Wise M., Nulsen P. E. J., David L. P., 
Sarazin C. L., Bautz M., Markevitch M., Vikhhnin A., 
Forman W. R., Jones C, Harris D. E., 2000, ApJ, 534, 
L135 

Mewe R., Kaastra J. S., Leidahl D. A., , 1995, Legacy, 6, 
16 

Pearson T. J., Readhead A. C. S., 1981, ApJ, 248, 61 
Pearson T. J., Readhead A. C. S., 1988, ApJ, 328, 114 
Peterson J. R., Fabian A. C, 2006, Phys. Rep., 427, 1 
Polatidis A. G., Wilkinson P. N., Xu W., Readhead 

A. C. S., Pearson T. J., Taylor G. B., Vermeulen R. C, 

1995, ApJS, 98, 1 
Rafferty D. A., McNamara B. R., Nulsen P. E. J., Wise 



M. W., 2006, ApJ, 652, 216 
Rebusco P., Churazov E., Bohringer H., Forman W., 2005, 

MNRAS, 359, 1041 
Russell H. R., Fabian A. C, Sanders J. S., Johnstone R. M., 
Blundell K. M., Brandt W. N., Crawford C. S., 2010, MN- 
RAS, 402, 1561 
Sanders J. S., 2006, MNRAS, 371, 829 
Sanders J. S., Fabian A. C, 2006a, MNRAS, 371, 1483 
Sanders J. S., Fabian A. C, 2006b, MNRAS, 371, L65 
Sanders J. S., Fabian A. C, 2007, MNRAS, 381, 1381 
Sanders J. S., Fabian A. C, Taylor G. B., 2005, MNRAS, 
356, 1022 

Sanders J. S., Fabian A. C, Taylor G. B., 2009, MNRAS, 
396, 1449 

Schneider D. P., Hall P. B., Richards G. T., Strauss M. A., 
Vanden Berk D. E., Anderson S. F., Brandt W. N., Fan 
X., Jester S., Gray J., Gunn J. E., SubbaRao M. U. e. a., 
2007, AJ, 134, 102 

Shepherd M. C, Pearson T. J., Taylor G. B., 1994, in Bul- 
letin of the American Astronomical Society Vol. 26 of Bul- 
letin of the American Astronomical Society, DIFMAP: an 
interactive program for synthesis imaging., pp 987-989 

Simionescu A., Werner N., Bohringer H., Kaastra J. S., 
Finoguenov A., Briiggen M., Nulsen P. E. J., 2009, A&A, 
493, 409 

Taylor G. B., Fabian A. C, Gentile G., Allen S. W., Craw- 
ford C, Sanders J. S., 2007, MNRAS, 382, 67 

Werner N., Simionescu A., Million E. T., Allen S. W., 
Nulsen P. E. J., von der Linden A., Hansen S. M., 
Bohringer H., Churazov E., Fabian A. C, Forman W. R., 
Jones C, Sanders J. S., Taylor G. B., 2010, MNRAS, 407, 
2063 

Werner N., Sun M., Bagchi J., Allen S. W., Taylor G. B., 

Sirothia S. K., Simionescu A., Million E. T., Jacob J., 

Donahue M., 2011, ArXiv e-prints 
Whyborn N. D., Browne I. W. A., Wilkinson P. N., Porcas 

R. W., Spinrad H., 1985, MNRAS, 214, 55 
Xu W., Readhead A. C. S., Pearson T. J., Polatidis A. G., 

Wilkinson P. N., 1995, ApJS, 99, 297 



© 2011 RAS, MNRAS 000, ??-[II] 



